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The magnetocrystalline anisotropy energies (MAEs) of the ferromagnetic metals bcc Fe, fee and 
hep Co, and fee Ni have been ealculated by using the ab initio tight-binding method. Disentangling 
the strong correlation among the d orbitals with the Hamiltonian in the loeal spin-density approxi- 
mation, we have investigated the orbital polarizations induced by the Hubbard U and Racah B. The 
experimental MAE of fee Ni is found with the value of U close to that determined from experiments 
and used in other theories. With the optimized values of U and J, both the MAEs and the orbital 
moments for Fe and Co are in close agreement with experiment. 
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Obtaining the magnetocrystalhne anisotropy energies 
(MAEs) of Fe, Co, and Ni from ah initio calculations 
within the local-spin-density approximation (LSDA) to 
densityfunctional theory is of considerable current inter- 
est HE II II mi- From various high-quahty LSDA cal- 
culations, the best case is Fe, where the computed values 
differ from experiment by a factor of about 2. The result 
for hep Co is far worse, and for Ni, the sign is not even 
correct. The effect of the so called spin-other-orbit cou- 
pling is far too small to bring theory and experiment into 
accord @- The discrepancy between theory and exper- 
iment, especially in the case of Ni, is usually attributed 
to the LSDA. 

The LSDA predicted MAEs and orbital moments can 
be improved for Fe and Co "5, 't'I by introducing the 
Brooks' orbital polarization (OPB) term j^il which mim- 
ics Hund's second rule. However for Ni, the predicted 
easy axis is still wrong Pj. In OPB, OP is driven by the 
Racah parameter B, with an energy functional related 
to the orbital moment (L) given by Ai?oPB = ~\B{L)'^ 
. It was argued that the key parameter responsible for 
the exchange-correlation enhancement of the orbital mo- 
ments in solids is the Hubbard U rather than the intra- 
atomic Hund's second rule coupling 0. Recently, the 
experimental MAEs of Fe and Ni have been obtained^ 
in the LDA-|-t/ method ITl with the noncoUinearity of 
intra-atomic magnetization included. However, the au- 
thors found the MAE of fee Ni to be a very rapidly 
varying function of U (from —50 to 60 /xeV/atom). A 
slight change of the value of U O.leV) may predict 
the wrong sign. Given this sensitivity, it is highly de- 
sirable to disentangle the intra-atomic strong correlation 
with the Hamiltonian in the LSDA and therefore to clar- 
ify the effect of OP in first principles calculations. This 
is the purpose of the present paper. 

The basic features of the electronic structure of Fe, Co, 
and Ni can be understood on the basis of their two types 
of valence-electron orbitals The extended s, p (,and 
/) orbitals should be well described by the LSDA. The 
fairly localized d orbitals, for which the electron-electron 
interaction is between the localized and itinerant limits. 



are not adequately dealt with in the LSDA 11]. In order 
to disentangle the strong correlation with the LSDA, we 
express the Hamihonian as H = TJLSDA + fjC _^ fyso 
Here TJLSDA jg standard LSDA Hamiltonian, H'^ is 
the correlation within the d orbital subspace and H^'~' 
the spin-orbital interaction. 

We start with the LSDA Hamiltonian in the orthogo- 
nal representation of the tight-binding linear muffin-tin 
orbital method in the atomic sphere approximation (TB- 
LMTO-ASA) 113, 
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where C, A, and 7 are the self-consistent standard po- 
tential parameters. Because the electron-electron inter- 
action is included in H'^ , the on-site diagonal matrix ele- 
ment of the d orbital is replaced with Cd = (C^ +C'j^)/2. 
S'^{k.) is the structure constant matrix in the orthogo- 
nal representation with k running over the Brillouin 
zone (BZ). ijLSDA ig block-diagonal in the spin index a 
along the magnetization direction. The spin-orbit cou- 
pling matrix elements for d orbitals are calculated in the 
last iteration of the self-consistent field procedure and 
H^'^ is treated in the usual single-site approximation [T^ . 

Similarly to the LBA+U method [lfl,li3, we treat the 
screened interaction among the intra-atomic d-orbitals in 
the Hatree-Fock approximation (HFA), 
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where J7, 

is the on-site d occupation matrix in the spin-orbital 
space. Ea is the average interaction without spin and or- 
bital polarization. Umim3m2nn are determined by three 
Slater integrals Fq, F2, and ^4^^, which are linked to 
three physical parameters: the on-site Coulomb repulsion 
U = Fq, exchange J = ji{F2 + F4), and Racah param- 
eter B = ^(9^2 — 5F4). In terms of U and J, Ea is 
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expressed as ^Un"^ — t/+4j where is the on-site 

d-orbital occupation. The ratio F4/F2 is. to a good accu- 
racy, a constant ~0.625 for d electrons [III, which leads 
to the estimation B w O.llJ. The interaction energy 
E'^'^, which is rotationally invariant with respect to the 
basis, leads to an effective potential i?^ acting on the 
d-orbital subspacc. 

The spin polarization (SP) in LSDA is generically close 
to the Stoner concept with an energy related to spin mag- 
netization m of AEgp^^{m) = —^I{ni)m'^, where /(m) 
is of the order of 1 eVj^l- In the HFA, the average spin 
splitting for d electrons is driven by / = ■^{U + 4 J) , with 
OP determined hy UcS = U - J In the hmit B = 0, 
Umim^m^mi Only involvcs two spherical harmonics 
and U = J = I is approximately equivalent to the LSDA 
[3. Even with Ucb = 0, there is no simple relation be- 
tween E'^'^ and the orbital moment practice, 
the problem involving the OPs induced by the Hubbard 
U and Racah B can be solved numerically by working di- 
rectly with the site-diagonal elements of the occupation 
matrix. 

The MAE is calculated by taking the difference of two 
total energies with different directions of magnetization 
(MAE = Em — Eqqi for cubic structures and MAE 
= -E-ioio ~ -£'0001 for hep structure). The total energies 
are obtained via fully self-consistent solutions of H, with 
the double counting corrections to the total energy in- 
cluded. For the k-space integration, we use the special 
point method ^3 '^ith a Gaussian broadening of 50 meV 
We use 100'^ sampling points in the BZ for cubic 
structure and 100 x 100 x 56 points for hep structure. 
We have also included the occupation number broad- 
ening correction terms to the ground-state total energy 
[23 ■ Numerical convergence has been tested against the 
number of k-points and Gaussian broadening. We first 
calculate the electronic structure self-consistently using 
the scalar relativistic TB-LMTO-ASA method. Then we 
construct the Hamiltonian H , using the spin-orbital cou- 
pling constants corresponding to the d band center l] . 

Considering the fact that the spin moments are well de- 
scribed by the LSDA, for a particular value of U, we have 
chosen the parameter J such that the magnetic moment 
maintains the theoretical value from the LSDA without 
spin-orbit coupling. The calculated spin moments are 
almost independent of B. Because the strong correla- 
tion U and J are entangled with the LSDA potential in 
the LDA-I-J7 method \Xm, the dependence of the mag- 
netic moment on U and J is not clear. Moreover, the 
energy E"'' defined in the LDA+C/ method pil llll| is 
not zero even without any SP and OP. This may render 
the interpretation of the delicate MAE dependence on U 
quite difficult. In their LDA-I-J7 calculations of Fe and 
Ni, Yang et al. scanned the (C/,J) parameter space 
and obtained the path of U and J values which hold con- 
stant the theoretical magnetic moment aligned along the 
(001) direction. J increases with U in their parameter 
path, in contrast to the basic concept ^^"^"^ ~ Id- The 
SP and OP are treated on the same footing in our HFA 



scheme. The calculated MAEs versus the values of U are 
depicted in Fig. ^ The corresponding orbital moments 
are presented in Fig. 12 
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FIG. 1: MAEs of bcc Fe, fee Co, hep Co, and fee Ni as a 
funetion of Hubbard U. Two eurves in spdf basis are plotted 
for each case, one with B — O.llJ, the other with B = 0. The 
experimental values are indicated by the horizontal dotted 
lines (1.4, 65, and —2.7 /leV/atom for bcc Fe, hep Co and 
fee NiI2, and 2.0 ^eV/atom for fee Co [U). The MAEs 
calculated in spd basis with B — O.llJ are also presented. 

The ASA does not significantly affect the accuracy 
of MAE. In fact, the differences in the MAEs and or- 
bital moments calculated by the LMTO-ASA method 
and those calculated by the full potential (FP) LMTO 
method are negligible if the same partial wave expan- 
sion Zinax = 2 is used. In cubic structures, the difference 
between the MAEs with /^ax = 3 and Zmax = 2 is very 
small (< 0.1 /ieV/atom)P|. However for hep Co, the 
MAE changed sign when angular moment /„iax increased 
from 2 to 3jl||. The MAE calculated in the spdf LMTOs 
is closer to the recent accurate result calculated from FP 
linearized augmented plane-wave (LAPW) method % 
Because the LMTOs in the spdf basis are more complete 
than the ones in the spd basis[l3|. we regard the results 
in spdf basis more reliable. 

In the limit B — and U J , our calculated MAEs 
and orbital moments (the left ends of the B = Q plots) 
lie in the range of the recent high quality LSDA values 
P, 12, 0,0, [5, 6J. In particular, the predicted signs of the 
MAEs for hep Co and fee Ni are wrong in the spdf basis. 
The calculated orbital moments are about 40% smaller 
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FIG. 2: Orbital magnetic moment for bcc Fe, fee Co, hep Co, 
and fee Ni along the experimental easy axis as a function of 
U . The plots are labelled the same as FIG.0 



than the experimental values for Fe and Co, while for 
Ni, it is very close to the experimental value. Since the 
LSDA is even poorer in Ni than in Fe and Co[l3|, the 
agreement must be accidental. 

By turning on the orbital polarization induced by the 
Racah _B at J7 « J, the calculated MAEs for Fe and Co 
are in much better agreement with the experimental data. 
This is similar to the LSDA+OPB calculations HQ. Par- 
ticularly, the predicted easy axis is correct for hep Co ■ 
The calculated MAE in s-pdj basis for hep Co is quite 
close to that of the LMTO-ASA calculation with OFBQ, 
while our spd MAE is quite close to the FP-LMTO+OPB 
resultQ. Interestingly, as shown in Fig. [21 the enhance- 
ment of orbital moment due to OP induced by Racah 
B is in excellent agreement with the OPB calculations 
despite the forms of E'^° and AE'qpb being quite dif- 
ferent. We suggest that the widely used OPB can be 
brought 'precisely into accord with the unrestricted HFA 
with A£;opB replaced hy E""^ at U ^ J in Ea re- 
placed with and n^). Only with the OP induced by 
B, the predicted sign of MAE for fee Ni is wrong. 

We now study the effect of OP induced by U with 
UeS > 0. Similarly to the OP induced by B, as shown in 
Fig. 121 the OP induced by U enhances the orbital mag- 
netic moments. For both Fe and Co, the MAEs change 
very smoothly and monotonically with increasing U. For 
Ni, the MAE is about zero when U « 1.3 eV. It de- 



creases with increasing U, then making a flat region with 
MAEs:! -1 eV from U = 1.7 to U = 2.5. After the flat 
region, the MAE decreases with increasing U, reaching 
the experimental value a,t U — 2.95 eV. The MAE pre- 
dicted here changes very smoothly with increasing U and 
with correct sign when [/ > 1.3 eV, without the strong 
sensitivity observed in the LDA+C/ calculations j^. 

TABLE I: U and J (in eV) corresponding to experimental 
magnetic anisotropy energy (in spdf basis with B — O.llJ). 
The calculated orbital moments {fj.B / atom) along easy axis 
are compared with the experimental data P]. 





bcc Fc 


fee Co 


hep Co 


fee Ni 


u 


1.15 


1.41 


1.77 


2.95 


J 


0.97 


0.83 


0.75 


0.28 


h 


0.087 


0.123 


0.150 


0.064 


expt. 


0.08 




0.14 


0.05 



When U > 2.5 eV, the two MAE curves of B = 
and B = O.llJ are almost indistinguishable for fee Ni. 
Thus we conclude that it is the Hubbard U that is fully 
responsible for bringing theory into accord with experi- 
ment. For Fe and Co, both OPs induced by U and B are 
needed to produce the experimental MAEs. As shown in 
Table m for bcc Fe and hep Co, the optimized U and J 
almost simultaneously give the experimental MAEs and 
orbital moments. The predicted orbital moment for fee 
Ni is slightly higher, but quite acceptable. The optimized 
values of U and J and their trend from Fe to Ni are very 
similar to those determined from experiments and used 
in other theories 0, The optimized values of U 

for fee Co and hep Co are close but not the same. This 
may be due to the fact that the experimental MAE of fee 
Co was extracted from measurement on supported films 
[2lj | or that the calculated MAE of hep Co is not fully 
converged even with /max = 3. 

It was conjectured that the failure of the LSDA to pre- 
dict the MAE of fee Ni is related to the band structure 
along TX direction 0. As shown in Figj^l similarly to 
the LSDA bands[ll, five bands cross the Fermi energy 
almost at the same k point when J k, JJ — 0.9 eV and 
5 = (The bands for B = O.llJ are quite similar). One 
of the d bands is just above the Fermi level at the X 
point and results in the appearance of small X2 pocket 
on the Fermi surface 0,0, which has not been found ex- 
perimentally. Increasing the valence electrons and thus 
pushing down the band corresponding to the X2 pocket, 
Daalderop et al. 1] found the correct easy axis for fee 
Ni. We find that the X2 pocket disappears at J7 « 1.5 
eV, corresponding to the start point of the correct sign of 
the MAE. With increasing [/, this band is further pushed 
down and when L/ ~ 3 eV, it is below the Fermi energy 
by about 0.2 eV. In our approach, the disappearance of 
the X2 pocket is a natural result of the OP induced by 
Hubbard U. 

As a delicate property, the MAE is naturally expected 
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FIG. 3: Band structure of fee Ni along Az = 27r/a(0, 0, 1) with 
the spin quantization axis in the (001) (solid lines) and the 
(111) (dashed lines) direetions. Here 0.4 </ < 1.0. I = 1.0 
eorresponds to the Xz point. 



to depend on the delicate changes of the band structures. 
We have applied the current scheme to the parametrized 
TB models of Fe and Ni fitting to the APW bands in 
the LSDA |2^. The trend and optimized values of U 
and J are found to be very similar to the ab initio TB 
calculations. This underlines the importance of how to 
treat the intra-atomic strong correlation. 

In summary, we have calculated the MAEs of Fe, Co, 
and Ni from the ab initio tight-binding total energies. 
Disentangling the strong correlation among the intra- 
atomic d orbitals with the Hamiltonian in the LSDA and 
therefore treating the SP and OP on the same footing, 
we have solved the long-standing notorious problem of 
the MAE of fee Ni. The discussions on the OPs induced 
by Hubbard U and Racah B can shed light on future first 
principles and model TB calculations. How to calculate 
the interaction parameter U in metallic environment di- 
rectly from first principles is still an open problem. 
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under the Competitive and Sustainable Growth Pro- 
gramme. 
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